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Abstract 

Local non-Hermitian potentials V{x) ^ V*{x) can, sometimes, generate stable bound 
states iplx) at real energies. Unfortunately, the idea [based on the use of a non-Dirac 
ad hoc metric Q{x, x') 7^ 6{x — x') in Hilbert space] cannot directly be transferred to 
scattering due to the related loss of the asymptotic observability of x [cf. H. F. Jones, 
Phys. Rev. D 78, 065032 (2008)]. We argue that for smeared (typically, non-local 
or momentum-dependent) potentials V V'^ this difficulty may be circumvented. 
A return to the usual (i.e., causal and unitary) quantum scattering scenario is then 
illustrated via an exactly solvable multiple-scattering example. In it, the anomalous 
loss of observability of the coordinate remains restricted to a small vicinity of the 
scattering centers. 
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1 Introduction 



An intuitive understanding of various physical aspects of scattering can be facilitated 
when one turns attention to simplified, one-dimensional schematic models of experi- 
mental setup. An extremely exciting toy-model scenario has recently been proposed 
and analyzed in ref . |1] . One of the most elementary and popular delta-function po- 
tentials Vo{x) = —a 6{x) has tentatively been combined with a remote non-Hermitian 
interaction. The purpose of this Gedankenexperiment has been formulated as a study 
of an "interface" between Hermitian and non-Hermitian components of the potentials 
exemplified by the superposition 

V{x) = Vo{x) + if3[6{x-L)-6{x + L)], L > 1. (1) 

Our present text offers an immediate continuation of this project. We feel motivated 
by the occurrence of many open questions in such a setting. In particular, the 
results of ref. [1] indicated that it might be rather difficult to keep a non-Hermitian 
interaction model short-ranged and compatible with the standard requirements of a 
local and causal physical interpretation of incoming and/or scattered waves. 

Our present answer to these compatibility questions will be predominantly af- 
firmative. More precisely, we shall emphasize that the number of problems which 
arose during the analysis of potential ([T]) may be attributed to its strict locality. We 
shall propose and advocate the replacement of the strictly local interaction operators 
V = V{x) by their slightly smeared descendants tractable as weakly momentum- 
dependent operators. For illustration purposes we shall use interactions V ^ V'^ 
given by eq. (IT^ in section [2] below. 

The latter choice of amended model will preserve its maximal similarity with 
the original potential of ref. p]. Firstly, the role played by the real "measure of 
non-Hermiticity" j3 in eq. ([1]) will be transferred to another real coupling constant g. 
Secondly, the distance L between the strictly localized interaction points of eq. ([T]) 
will be replaced by a variable integer A/" representing a separation distance between 
two not entirely local, "smeared" domains of support of our interaction V. 

We decided to parallel the majority of quantitative results of ref. [I| by their 
close and explicit analogues. The major problems resulted from the manifest non- 
Hermiticity of the interaction which implies, for eq. ([T]) at least, the necessity of a 
drastic change of the concept of the coordinate. In ref. [2] this mathematical result 
has been identified as a source of a deep conflict between the use of x in eq. ([1]) (i.e.. 
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in the "input" definition of tfie interaction) and, simultaneously, in tlie asymptotic 
boundary conditions for tlie one-dimensional scattering. 



It is necessary to keep in mind that the core of this conflict does not lie in the 
formalism of quantum mechanics itself. Formally, no problems occur since all the 
unitary transformations of a given model (cf., e.g., a non-local free- motion example 
given in section 5 of ref. [3J) must lead to equivalent physical predictions. 

The differences in predictions can only occur when non-equivalent definitions of 
the dynamics are being compared. This is precisely in this sense that the difficulties 
emerged in refs. [H [2] where a simultaneous validity of both the definition ([T]) of the 
local physical interaction V V"^ and an a priori assignment ([2]) of the usual physical 
meaning to the local free waves exp ±iKx in asymptotic domain has been required. 

A key to the resolution of this misunderstanding has been described in our paper 
|1] where we showed that there exist non-Hermitian models where the simultaneous 
validity of the phenomenological postulates ([1]) and ([2]) can be achieved after a certain 
modification of their respective forms. In this sense, our present paper will just amend 
and strengthen the argumentation of ref. [1]. Indeed, in a way emphasized by the 
note added in proof in [T], our old model "still involved a departure from standard 
quantum mechanics at large distances" . 

The present final resolution of the conflict between the locality of forces and 
waves will rely on a nontrivial extension of the class of interactions accompanied by 
an enhancement of efficiency of necessary mathematics. These technical details will 
be described in sections [3] and H] and in two Appendices. In nuce, the Runge-Kutta 
coordinate-discretization method [5j will be shown superior, for the given purpose, 
to the usual perturbation expansions as employed, e.g., in ref. [1]. In the latter study 
of eq. ([T]), for this reason, both the inverse length 1/L and the variable coupling 
constant (3 had to be assumed small. In contrast, the variability of both our present 
parameters g and Af will be, within their respective physical ranges, unrestricted. 

The presentation of our explicit scattering solutions in section H] will confirm the 
full consistency and unitarity of the scattering in our amended class of non-Hermitian 
models. In summary (section several comments will finally be added clarifying 
the proposed changes of theoretical perspective in a broader, less model-dependent 
context. 




x > 1 . 



(2) 
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2 Towards the short-ranged non-Hermiticities 



2.1 Runge-Kutta discretization 



— 3-^ + V'(x) = , xG(— 00,00) (3) 



We may treat any one-dimensional Schrodinger equation 

with a local and real or complex potential V{x) as a continuous /i ^ limit of 
its difference-equation approximation defined along the Runge-Kutta doubly infinite 
lattice of discrete coordinates x = = k h, k = 0, ±1, . . ., 

Approximate wave functions may be then constructed via the methods of linear 
algebra re-parametrizing, incidentally, the real energies E = {2 — 2cosip)/h'^ in terms 
of a real angle (f = (f{E) G (0, vr). The scattering boundary conditions ([2]) may and 
should be re-written in their discrete version, 

^ ^ ( e'"''^ + Re-'""^ , m<-M<-l, 

-ipiXm) = < . (5) 

[ Te*'"'^, m > M - 1 . 

One could easily discretize the ultralocal non-Hermitian toy model ([1]) and confirm 
the discouraging conclusions, formulated in ref. jjj, that one can "no longer talk 
in terms of reflection and transmission coefficients" so that "the only satisfactory 
resolution [of dilemmas] is to treat the non-Hermitian scattering potential as an 
effective one, and work in the standard framework of quantum mechanics, accepting 
that this effective potential may well involve the loss of unitarity" [T]. 

The loss of unitarity need not necessarily be perceived as a weakness of the theory, 
especially when one deals "with a subsystem of a larger system whose physics has not 
been taken fully into account" [Ij . In this sense one may perceive eq. ([1]) with a local 
non-Hermitian interaction V{x) 7^ V*{x) as an "effective theory", i.e., as an incom- 
plete picture of physical reality. This philosophy finds interesting phenomenological 
applications ranging from classical optics [B] or models with supersymmetry [7] up 
to the manifestly non-unitary scattering models in quantum phenomenology [8] and 
up to the descriptions of open systems in nuclear and solid state physics [9] and in 
quantum cosmology [TU]. Nevertheless, in an alternative, theoretically much more 
ambitious approach to the localized non-Hermitian scattering potentials one should 
insist on the conservation of a suitable current and, hence, on the strict unitarity of 
the scattering realized by the asymptotically observable free plane waves. 
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2.2 Nonlocal updates of potentials 



The first (and the least sophisticated) quantitative non-Hermitian model satisfy- 
ing the above requirements has been constructed in our paper [4j. We replaced 
equation (jl]) by its generalization where the interaction operator V acquired the 
nearest-neighbor form, 



+ 



+ Vk,k+i i^{xk+i) + Vfc,fc ipixk) + Vfe,fc-i ^(a;fc-i) = Etp{xk) ■ (6) 

After a re-scaling of the Hamiltonian H = —dp/dx^ + F by an inessential numerical 
factor we obtained 



H = -^ + V, -A 



2 -1 
-1 2 -1 

-1 2 -1 
-1 2 



(7) 



and choose the following, minimally non-local potential 



— c 

c -h 

h 


—a 


a 


-h 

b -c 
c 



The resulting multiparametric Hamiltonian H = —A + V^"'''^'^'"'^ ^ remains 
manifestly non-Hermitian in the "friendly" Hilbert space H^^^ endowed with the 
usual inner product 

in n^^K 



(9) 



Note that the summation would only be replaced by the integration in the continuous 
limit h ^ 0. Now, the key point is that the same operator H may be found Hermitian 
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after one moves into another Hilbert space Ti.^^^ . In the latter space the definition of 
the inner product must be different and more general, 

{m^'^ = Y.Y.r{xk)ek,nHxn) = {m\<i>)--={{m m n^'^. m 

k n 

The "non-Dirac metric" matrix = 0^ must only remain positive definite and 
compatible with the Hamiltonian in question [TT], 

H^Q = QH. (11) 

In the notation of ref. [12] one writes H = and speaks about a "quasi-Hermiticity" 
|llj or "pseudo-Hermiticity" [13] or "cryptohermiticity" [H] of the Hamiltonian. In 
this context the core of the message delivered by our paper [1] was that there exists 
a metric-operator matrix remains compatible with our interaction 

model ([HD as well as with the asymptotic observability of Runge-Kutta coordinates 
Xk- This matrix has the following compact and fully diagonal form, 



0-5 

0-3 






03 

05 



(12) 



Its elements are given by closed formulae, 

e±i = (1 ± a)(l - b^){l - c^){l -d^)... , 

^±3 = (l±a)(l±6)2(l-c2)(l-rf2)... , 
^±5 = (1 ± a)(l ± 6)2(1 ± c)2(l -d^)... . 

One arrives at a causality-observing physical picture of scattering based on a clear 
separation of the "in" and "out" solutions not only in Hilbert space Ti.^^^ but also in 
Hilbert space 

In our subsequent paper [15] the next step has been made. In the spirt of eq. ([1]) 
we simulated the existence of several separate point interactions. Unfortunately, the 
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construction of the metric only remained feasible under a very specific, left-right 
symmetric arrangement of the set of interaction centers. Sometimes, this type of 
symmetry is being called PT— symmetry, for reasons and with motivations which 
are thoroughly explained elsewhere [TB] . 

Our present continuation of development of the multiple-scattering idea will be 
based on a return to asymmetric models, allowing an independence of arrangement 
of several spatially separated scatterers. Paradoxically, the transition to asymmetric 
realizations of the set of interaction centers will be accompanied by a simplification 
of analysis of their mutual interference. 

2.3 Limiting transition to continuous coordinates h ^ 

For a quantitative specification of the extent of nonlocality induced by multiparamet- 
ric matrices V^"'''''^'"'^ of eq. ([8]) let us start from the simplest, coordinate-independent 
model where a ^ 6 ^ c ~ . . .. Then, the limiting transition to h = con- 
verts operator into the first power of the momentum, \/('^'"'"'---) ~ d/dx. 
In the subsequent step one may re-introduce a weak coordinate-dependence (with 
a ^ b ^ . . .) and evaluate the continuous limit perturbatively. Locally, the limit 
h will preserve the same leading-order approximate proportionality of the 
coordinate-dependent potential to the momentum. 

Admitting an unconstrained variability of the parameters in matrices V^"'''''^'"'^ we 
obtain some less trivial coordinate- and momentum-dependent operators. For the 
sake of brevity let us restrict similar considerations solely to the models with just 
a few nonvanishing coupling parameters. Then, the limiting transition h —>■ will 
certainly lead to point interactions. Their explicit definition will be given precisely 
by the matching of the wave functions. Just a slightly more complicated alternative 
to the delta-function point-interaction model ([T]) of ref. will be obtained. Our 
Appendices A and B may be consulted for illustration of some technical aspects of 
such a type of matching recipe. 

Our illustrative toy potential ([8]) has not been too well designed for phenomeno- 
logical purposes since it did not allow us to remove the spatial asymmetry from the 
related metric matrix f[T^ . 

0-k _ {i-a){i-bni-cr... 

Ok (1 + a)(l + 6)2(1 + c)2... ■ ^ ' 

The effect of the localized non-Hermiticity in H remained long-ranged. 
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3 Toy model 



Equation ( fT3l) indicates that the flow of the probabihty is different to the left and to 
the right of the scattering center. A weaker form of this shortcoming characterizes 
also the PT— symmetric models of ref. [15j where the metric remained rescaled (i.e., 
non-Dirac, 0^^^ 7^ 1) along the spatial interval(s) separating the individual scat- 
terers. This encouraged us to perform a series of computer-assisted trial-and-error 
experiments leading, at the end, to our present interaction-matrix candidate 
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large gap. 


-9 
9 9 
-9 









27V + 1 
columns 








where each scatterer is simulated by three-dimensional submatrix. Although our 
particular model (|T4l) comprises just two localized interaction centers at a;±(Ar+2), 
we shall not consider three (like in eq. ([1])) or more individual scatterers because 
such a generalization would remain routine, not necessitating any significant further 
technical improvements of our method. 

3.1 Metric Q^^X) with localized anomalies 

Using heuristic arguments we arrived at ansatz and studied the scattering solu- 
tions. Schrodinger equation with the smallest gaps M has been studied first of all. A 
sample of these calculations may be found collected in Appendices A and B below. 
They demonstrate that one of specific merits of eq. (|T^ lies in a maximal simplicity 
of necessary algebraic manipulations. 
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The second merit of the choice of eq. f|T^ can be seen in its generic character. One 
can add several further interaction submatrices of the same form without worsening 
the feasibihty of the calculations. On this background, without any real loss of 
generality we restricted our attention just to the first nontrivial example which is 
characterized by the occurrence of the mere two remote centers of interaction. 

We decided to construct all the eligible metric matrices as linear-algebraic so- 
lutions of eq. (fTTj) . After we imposed the condition of the compatibility of with 
the asymptotic observability of the coordinate we revealed that our present models 
can be assigned the diagonal metric operators of the same doubly infinite 
diagonal matrix form 



1 












1 

1 
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1 
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= 2M+l 
("distance") 




1 



This metric differs from the Dirac's = / solely at the centers of the non- 

vanishing three-by-three submatrices simulating the non-Hermitian point-like scat- 
terers. 

3.2 Single-center limit /i^OatAA=— 1 

The picture of the scattering as offered by our toy potential (fT^ and by the related 
metric matrix fllSp depends on the Runge Kutta discretization length /i > 0. Once 
we demand that a measured distance between two scattering centers is a macroscopic 
constant L, our parameter M must grow with the decrease of h as L/h. Vice versa, 
the use of an /i— independent M will only lead to a single-centered scatterer. In the 
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latter scenario the scattering is realized by a "quasi-local" potential. Its explicit 
specification will depend on the /;.— dependence of A/" = M{h). It remains compatible 
with L = whenever hM{h) —>■ for h ^ 0. Such a flexibility may make the 
interactions better suited for fine-tuning, say, of the strength of non- localities and/or 
of the extent of the violation of conservation laws at short distances, etc. 

For illustration purposes let us pick up the elementary example of Appendix A. 
Relaxing the specification of some concrete asymptotic boundary conditions let us 
re-interpret its "distant" wavefunction components as an arbitrary free wave ip{x). 
At X < X-2 or X > X2 this yields the coincidence of symbols 

U-ni = IpiX-rn) ■= V'-m'"''' , Lm = 1p{Xm) ■= Ip+m"^ , m > A/" + 3 = 2 

respectively. Next, the first and last matching condition extend both the latter 
assignments by one more step, 

U.i = + Li = {1 + g) ij{x+i) . 

Finally, with ipo = ipi^o) we arrive at the three dynamically nontrivial requirements 



2cOS(y9 —1 




f/-l 






-1 + 5(2 2cosv2 -l+g^ 











—1 2 cosy) 




^1 







which define our wave function, implicitly, near the origin. Tentatively, we may 
Taylor-expand 

^i^f ^) 2h^' + 2h^ij" + . . . , (1 + g)-^U^i = ij-hij' + h^i)"/2 + . . . , 

7/'o = ^, {l + g)-^U = ij + hij' + h^ij"/2 + ... , ^5^"^'^ =^/' + 2/l^' + 2/^V" + ■■■ 
and insert these approximants in eq. f|T6l) . yielding 

-{l-g){'4)- 2hil}' + 2h'^^") + 2 cos(^ - hip' + h^i)" /2) - (1 - 5) ^ = 0{h^) , 

-(1 + 5) - ht^' + h^iP"/2) + 2cos(/?V^ - (1 + 5) (V^ + hi)' + /iW2) = 0{h^) , 

-{l-g){i) + 2hi>' + 2/iV') + 2cos^{i) + hip' + h^ip" /2) - (1 - 5) ^ = 0{h^) . 

According to these relations, Schrodinger equation V — E = ip" /ip would make 
quantity V large and positive when extracted from the combination of the first and 
third equation, or large and negative when extracted from the middle equation. This 
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means that our tentative assumption about the smoothness of wave functions near 
the origin leads to mathematical contradictions and must be abandoned. 

Let us now modify our assumptions, distinguish between the left and right wave 
functions and set A{x) = ip^x — h) and B{x) = il){x + h), i.e., 

(1 + g)-^U-i = A(0) := A , (1 + g)-^U = 5(0) := B . 

Naturally, 

while quantity ipQ acquires the two alternative first-order representations, 

ipo^ A + hA' ^ B -hB' . 

In the limit h ^ the latter relation yields 

A = B, A' = -B' . 

The insertion of our amended ansatzs in eq. f|T6l) leads just to the three alternative 
versions of the requirement of smallness of A = 0{h?), B = 0{h^) as well as of 
ipQ = 0{h). Thus, in the continuous-coordinate extreme our simplest Af = —1 
example degenerates to the opaque- wall-barrier dynamics generated by an additional 
Dirichlet boundary condition V^(0) =0. 

We see that the role of non-Hermiticity is, in our model with Af = —1 at least, 
truly non-perturbative and dynamically highly influential. This conclusion may inde- 
pendently be confirmed by the inspection of the Af = —1 reflection and transmission 
coefficients given in Appendix A. We believe that also beyond this concrete example, 
at least some of its features will survive a transition to more-center models and/or 
to the two-center models at large separation distances Af{h) = 0{l/h). 

4 The unitarity of the scattering at any J\f 

The first encouraging surprise encountered during the inspection of the discretized 
metric (ITB]) is that it remains asymptotically diagonal in the coordinate representa- 
tion. This means that the asymptotic coordinate x remains observable. Moreover, 
the range of influence of individual non-Hermitian scatterers is shortened via eqs. f|T^ 
and f[T5|) . Thus, the only missing component of the whole picture are formulae for 
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the reflection and transmission coefficients, the determination of which may start 
from the hnear Schrodinger equation for discretized wave functions = iplxk) = ipk, 

Hi} = E^. (17) 

In its hght the vahdity of boundary conditions ^ can be prolonged to all the sub- 
asymptotic free-motion domain, 

In parallel, for larger integers A/" = M^h) we may profit from adding another free- 
motion ansatz at the smaller subscripts, 

^jjk = C e'^'P + D e-'^"^ , \k\<J\f. (19) 

One should add that the study of the large distances A/" ^ 1 might be well motivated 
by its potential relevance in physics. In particular, its feasibility could offer a guide 
for simulation of macroscopic non-localities, the presence of which could, in its turn, 
lead to the violation of causality at small distances. In parallel it is important 
that the effect of our non-Hermitian V can be kept localized. This means that in 
contrast to virtually all of the published older models the simplicity of interaction 
(fT^ enables us to return to the "old-fashioned" definitions of the reflection coefficient 
R and transmission coefficient T. 



m 



> Ar + 3, 



4.1 The elimination of J\f from matching conditions. 

The second surprise offered by our example is that the matching remains easy even 
for remote interactions with A/" ^ 1. In order to show this, let us now assume that 
the distance 2A/'+ 1 between two three-dimensional interaction submatrices in (IT^ 
is arbitrary. We may abbreviate, in partitioned notation. 
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where denotes a null matrix (of dimension 2M + 1 ) and where are null column 
vectors. The superscripts denote transpositions (i.e., row real vectors). In such a 
notation one has to consider the following 2M + 7 matching conditions 
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and where = (—1, 0, . . . , 0) and = (0, . . . , 0, —1) are two (2A/'+ 1)— dimensional 
auxiliary row vectors. The other auxiliary "free-motion" submatrix F^^^ is tridiag- 
onal and (2A/' + 1)— dimensional. Its elements 2 cosy? along the main diagonal are 
complemented by the elements —1 which lie along its two neighboring diagonals. 



4.2 Exact solvability. 

What remains for us to demonstrate is that our model conserves the global or asymp- 
totic flow of probability, i.e., that one obtains |-Rp + |Tp = 1 in spite of the manifest 
non-Hermiticity of the Hamiltonian H. In this setting the final surprise comes with 
the observation that the reflection and transmission coefficients are obtainable in 
closed form. Even when the "distance" parameter Af is arbitrarily large, the use of 
ansatz (fT9l) reduces the original set of 2A/' + 7 matching conditions to the following 
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two independent matching conditions consisting of four items each, 



2 cos(p —1 — 


9 








-1 + g 2 cos 


i -I 


-1+9 





-1 - 


9 


2 cos V? 


-1 







-1 


2 cos 99 





U-M-3 










-M-2 + X-2 











-N~i + X-1 




















2cos</7 


-1 








-1 


2 cosy? 


-1-9 








-1 + 9 


2 cos (/9 


-l+g 










2 cos if 



Lm+i + Xi 
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Lj\f+4, _ 



Out of this octuplet of equations, the first and last hues can be solved, 

(1 + 9)x-2 = -gU-M-2 , (1 + g)x2 = -gLM+2 ■ 
This leads to the following two triplets of conditions 



2 cos if 
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-l+g'^ 2cos(^ 
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_ {l-g')Lj^+s 



Using the first and last equation we eliminate 

(1 - 9^) X-i = 9^ t^-AT-i + 9^U_^_3 , (1 - g^) Xl = 9^ Lj^+i + 9^Lj^+3 
The net result of these manipulations are the four relations 



2cos</7 —1 
—1 2cos</7 



U-M-l + 9^U^M- 

(1 - 



{l-g'')U^Af^2 
(1 - g'^)^-Af+i 



2cOS(/9 —1 

—1 2cos(^ 



-f'A/'+l + 9^Lj^+3 



{l-g')^lJM.^ 
(1 - g'^)L^+2 



which can be simplified to read 

(1 - g'') V'-AT = + 2g^U-M-2 + g^U^M-i , 
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(1 - 9'^) V'A/- = -t'Ar + 2g'^Lj^+2 + 9^Lx+4 , 

(1 - 5'^) i>j\r+i = Lj^+i + 9'^Lj^+3 ■ 

These equations represent the two alternative definitions of the sum C + D and of 
the difference C — D of the two unknown coefficients in ■0^, 

2{l-g^){C + D) cosAfip = A* (ip) + A{ip) (R + T), 

2 (1 - /)(C + D) cos(AA + l)v9 = B*{^) + B{^) {R + T) 
-2 i (1 - /)(C - D) sinAfip = A*{ip) + A{ip){R - T) , 
-2 i (1 - /)(C - D) sin(AA + l)ip = B*{ip) + B{(p) {R - T) 
where we abbreviated 

In the next step we ehminate C and D and express 



' sin(A/'+l)(^ sinjV(^ ' 



v{^p) ' C0s(jV+l)</7 COSj\f(p 

The required amphtudes and T are now found, in closed form, as the respective 
sum and difference of the latter two expressions. In the final step their probability 
conservation property 

|i?|2 + |r|2^i 

is easily seen. 

5 Summary 

Our main technical result is that via a discretization of the real axis of coordi- 
nates X (and using the matching method) an exact linear- algebraic solvability of our 
present model of scattering has been achieved. Constructively, the necessary unitar- 
ity requirement has been satisfied at the same time. Our model [containing several 
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spatially separated and strictly localized interactions which appear non-Hermitian in 
= L2{1R)] is being assigned the more or less unique Hilbert space of states 

^ ^{physical) ^J^g^g 

use of an anomalous inner product makes the Hamilto- 

nians (crypto) Hermitian. 

It is worth noticing that the metric operator which defines the inner product 
in '}{iphysicai) Yn.eie\Y differs from the usual Dirac's delta function locally^ viz., in a 
close vicinity of interaction points. This implies that the physical operator of the 
coordinate remains unmodified almost everywhere. An entirely consistent physical 
picture of scattering from multiple scatterers is obtained in this way. In contrast 
to the older models using non-Hermitian but strictly local potentials V"(x), the free 
motion between our present, slightly nonlocal individual non-Hermitian point-like 
scatterers remains undistorted. 

In conclusion let us re-emphasize that the motivation and inspiration of our 
present study of a simplified model of multiple scattering resulted from several 
sources. One of the most important ones has to be seen in the recent enormous 
growth of interest in the models of quantum dynamics of bound states which look 
manifestly non-Hermitian in L2{M) and/or in similar mathematical representations 
of the Hilbert space of states [16]. 

The key to success can be seen in the discovery of feasibility of a strictly physics- 
motivated transition to correct Hilbert space 

^{physical) U2]. Our present paper 
can be read as an implementation and advertisement of such an approach where one 
chooses a slightly more complicated input physics (i.e., in our case, a slightly nonlocal 
Hamiltonian H = T+V) and where one is rewarded by a perceivable simplification of 
mathematics. In particular, we saw that the resulting metric B in '}{(p^y'^^<^°-'-) differed 
from the unit operator just in a finite number of matrix elements in our model. 

We have only to repeat that our second, less abstract motivation grew from 
the emergence of several very recent studies of manifestly non-Hermitian models of 
quantum scattering [S]. Several issues may be addressed in this context. For example, 
in the less ambitious, effective-theory versions of these models (where one does not 
insist on the conserved probability) one can easily stay in the single, effective-theory 
Hilbert space Ti.^^^ Moreover, various additional dynamical assumptions (like the 
strict locality of potentials V{x)) may easily be incorporated in the similar pragmatic 
applications of the theory. 

In contrast, in the "fundamental" and unitary quantum theory a real challenge is 
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to be seen in the existence of a correlation between non-Hermiticity of a local V and 
the long-range non- locality emerging in [T7] . One can notice that this relationship 
seems highly model-dependent. In this sense, our present message can be read as 
a methodical encouragement. Basically, we found that whenever one broadens the 
class of the eligible potentials the latter model-dependence can be re-interpreted as 
an advantage. 

It should be remembered that the increase of the non-Hermiticity of H need not 
necessarily be correlated with the growth of non-localities in B obscuring the clear 
physical picture of scattering. We succeeded here in showing that both the non- 
localities occurring in V and O can be kept under control simultaneously. After all, 
one may note that in our present one-dimensional model with x & M the anomalies 
disappear "almost everywhere" in the continuous limit h ^ 0. 

In this manner our present text brought a rather surprising resolution of the puz- 
zle formulated in ref. [15] where we did not manage to get rid of the non-locality 
in a non-Hermitian model comprising several spatially separate scattering centers. 
Here we revealed that sometimes it makes good sense to sacrifice some inessen- 
tial symmetries of the model in order to preserve either its exact solvability or its 
phenomenological flexibility. We should note that the feasibility of our (computer- 
assisted) algebraic manipulations survived even the transition to unusually compli- 
cated point-interaction simulations by three-by-three matrices. 

In the context of physics good news concern, first of all, the possibility of an 
explicit construction of an optimal metric B in the physical Hilbert space. Its "op- 
timality" reflects the fact that with an obvious exception of the closest vicinities of 
the point-like interaction centers of our model, the metric B itself has successfully 
been forced to commute with the operator of the coordinate almost everywhere. This 
means that in contrast to intuitive expectations (supported even by some solvable 
models), the concept of coordinate and of an asymptotically free (i.e., measurable) 
motion of a quantized object can survive the emergence of a finite number of point- 
like non-Hermitian obstacles positioned arbitrarily along the real line. 

The latter observation allows us to declare that our model represents an illustra- 
tive example of a standard quantum system where the non-Hermiticity as well as the 
resulting non-localities (in both the metric B and in wave functions) remain confined 
to a very small part of the domain of the coordinates. This means not only that up to 
the singular points the coordinates remain measurable but also that the clear physi- 
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cal picture and consistent probabilistic interpretation of the non-Hermitian systems 
is naturally being extended to the multiple-scattering scenario. 
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Appendix A. Construction of the amplitudes for 
model (ED at A/'= -1 

In the special case of our toy model H^^"^^ at Af = —1 the analysis of the respective 
transition and reflection amplitudes T and R can be based on the explicit solution of 
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Schrodiner equation which degenerates, in an obvious manner and under the notation 
conventions of section |H to the following set of the five linear relations representing 
matching conditions near the origin, 
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Their solution may start from the first and last line giving 
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This enables us to consider just the three modified matching conditions 
—1 + g'^ 2 cos (p 
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—1 2cosv9 
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The first and last rows read 

(1 - g')iJo = Uo + g'U^2 =1 + 9' e"'^^ + {1 + 9' e^'^ R 

(1 - /)^o = Lo + g'L2 = {l+ / e^'n T 
so that their combination 

1 + = (1 + e"'^) (T - R) 

defines the difference between our two amplitudes as a complex number with unit 
norm, 

T-R = - — = e'" , A 



g"^ sin 2ip 



1 + iA ' 1 + fi'^ cos 2(p 

The remaining central matching condition can be given the form of an equation for 

the sum S = T + _R of the amplitudes, with the solution equal to another complex 

number with unit norm, 



T + R 



1/i 



1/i 



/i 
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[1-g^) sm2(p 



1 — 3g'^ — cos 2ip — g'^ cos 2ip 



This gives the two final formulae 



with the two respective properties 

4 = - e^") [e-'" - e'"') = 2 - e^^""^) - e'^'^"") 
which imply that 

|T|2 = 1. 

This means that in contrast to the observations made in some other non-Hermitian 
models [H [21 H] , the flow of probability is conserved so that the standard physical 
picture of the scattering does not require any modiflcations. 



Appendix B. Construction of the amplitudes for 
model ([HD at A/" = 



In place of the flve-dimensional matching condition of our preceding appendix let 
us now turn our attention to the family of nontrivial models where the two three- 
dimensional elementary-interaction submatrices are separated by a free-motion inter- 
val of the length 2M + 1. In the flrst nontrivial model with A/" = the nonvanishing 
submatrix of our interaction matrix is seven- dimensional, 
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In such a case one has to consider seven matching conditions of the form 
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The separate subset of the first and last matching condition is solvable as follows, 

il + 9)X-2^ -9U-2, {1 + g)x2^ -9^2- 

The backward insertion of these formulae leads to the quintuplet of the reduced 
matching conditions 
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Its first and last line define the other two correction components, 

(1 - x-i = / + C/-3) , (1 - XI = / (i^i + i^-a) 
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so that we are left with the three matching conditions 
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Their first and last item define the same quantity in two ways, 

{l-g^)ijo^Uo + g^ (C/-2 + 2 cos U^s) = C/q + / (2 C/_2 + C/-4) 

(1 - g^) i/jo = Lo + g^ (L2 + 2 cos <^ L3) = Lo + g^ (2 L2 + L4) . 
In effect, one can eliminate "00, 

(T -R) [1 + (2 e^^^ + e^^^ )] = [l + (2 e'^^^ + e"^^^ )" 

and specify the difference between T and i?, 

5(^(2 sin2(/9 + sin4(/9) 



1 + iX' 



A'- 



1 + 51^(2 cos 2(^ + cos iifi) 



Next, in a complete parallel to the previous construction, the sum E of T and R 
may and should be extracted again from the last and symmetrized middle item of 
our matching conditions, 

2C/_i + 2Li + 2/([/_3 + iv3) = Uo + Lo + g'' (2 C/_2 + 2 L2 + C/-4 + ^4) ■ 

After appropriate insertions this gives the similar formula as above, 

—2 sin ip + g'^{2 sin 2ip — 2 sin 2>(p + sin4(^) 



T + R^-^-^.^e^^' 



1 + ifj,' ' [1 — 2cos(/? + g'^(2cos2(^ — 2cos3(^ + cos4(/?)] 

The same argumentation as above confirms the validity of the identity 



i.e., of the same probability conservation law as above. 
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